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Abstract 

The deficit observed in the Galhum radioactive source experiments is interpreted 
as a possible indication of the disappearance of electron neutrinos. In the effective 
framework of two-neutrino mixing we obtain sin 2d > 0.03 and Am^ > 0.1 eV^. 
The compatibility of this result with the data of the Bugey and Chooz reactor 
short-baseline antineutrino disappearance experiments is studied. It is found that 
the Bugey data present a hint of neutrino oscillations with 0.02 < sin^ I'd < 0.08 
and Am? ~ 1.8 eV^, which is compatible with the Gallium allowed region of the 
mixing parameters. This hint persists in the combined analyses of Bugey and Chooz 
data, of Gallium and Bugey data, and of Gallium, Bugey, and Chooz data. 
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1 Introduction 



The observation of solar and very-long-baseline reactor neutrino oscillations due to the 
squared-mass difference ArngQ^ = (7.59 ± 0.21) x 10~^eV^ pLj and the observation 
of atmospheric and long-baseline accelerator neutrino oscillations due to the squared- 
mass difference Am^xM = 2.741q;25 x 10~^eV^ [2] give very robust evidence of three- 
neutrino mixing (for reviews of the theory and phenomenology of neutrino mixing, see 
Refs. [3|, m [5|, [6l [71 [HI m [TOl HI])- There are, however, some anomalies in the data of 
neutrino experiments which could be interpreted as indications of exotic neutrino physics 
beyond three- neutrino mixing: the LSND anomaly [12], the Gallium radioactive source 
experiments anomaly [13], and the MiniBooNE low-energy anomaly [H]. In this pa- 
per we consider the anomaly observed in the Gallium radioactive source experiments 
[T5l [T6| \T7\ [T8l [13], in which the Gallium solar neutrino detectors GALLEX [H] and 
SAGE [20] were tested by measuring the electron neutrino flux produced by intense ar- 
tificial radioactive sources placed inside the detectors. The Gallium radioactive source 
experiments measured a number of events smaller than expected. This deficit can be 
interpretecS as an indication of the disappearance of electron neutrinos due to neutrino 
oscillations [221 [221 1211 ■ Under this hypothesis, we analyze the data of the Gallium ra- 
dioactive source experiments in the effective framework of two-neutrino mixing, which 
describes neutrino oscillations due to a Am^ that is much larger than the solar and at- 
mospheric ones (see Refs. [5l[6l|TT]). We also study the compatibility of this interpretation 
of the Gallium radioactive source experiments anomaly with the data of the Bugey [25] 
and Chooz [26] reactor short-baseline antineutrino disappearance experiments. 



2 Gallium 

The GALLEX pjj and SAGE [20j solar neutrino detectors (see Refs. [HI H |5l [Q [3 [HI 
[9l[T0l[TT]) have been tested in so-called "Gallium radioactive source experiments" which 
consist in the detection of electron neutrinos produced by intense artificial ^^Cr and ^''Ar 
radioactive sources placed inside the detectors. 

The radioactive nuclei ^^Cr and ^''Ar decay through electron capture {e" + ^^Cr — >■ 
^^V -|- Ue and e~ + ^^Ar — > ^^Cl + z/g) emitting Ue lines with the energies and branching 
ratios listed in Tab. [H These neutrinos were detected through the same reaction used for 
the detection of solar neutrinos [7F] : 

i/e + ^^Ga-^ ^^Ge + e", (1) 

which has the low neutrino energy threshold E^(^^Ga) = 0.233 MeV. The cross sections 
of the z/g lines emitted in ^^Cr and ^^Ar decay interpolated from Tab. II of Ref. [2H] are 
listed in Tab. [H 

The ratios R of measured and predicted ''^Ge production rates in the two GALLEX 
^ Cr radioactive source experiment^ Crl [15] and Cr2 [I6], and the SAGE ^^Cr [13 [TH] 

^ Another possible explanation is that the theoretical cross section of the Gallium detection process 
has been overestimated [HI [5T] . 

^ As explained in Ref. [13], the values of R in Tab. [2] for the two GALLEX ^^Cr radioactive source 
experiments are different from those published in Refs. [HI [16], because of an improved reanalysis of the 



2 





51 Cr 


37 Ar 


[-K-t; V J 


747 


7'^') ZL97 




oil O 1 o 


B.R. 


0.8163 


0.0849 0.0895 


0.0093 


0.902 0.098 


a [10-46 p^2] 


60.8 


61.5 26.7 


27.1 


70.1 70.3 



Table 1: Energies (-Ejy), branching ratios (B.R.) and Gallium cross sections (a) of the 
z/e lines emitted in ^^Cr and ^^Ar decay through electron capture. The cross sections are 
interpolated from Tab. II of Ref. [2H]- 





GALLEX 


SAGE 




Crl Cr2 


51 37^1- 


R 




1.00 ±0.10 0.81 ±0.10 


0.95 ±0.12 0.79 ±0.10 


radius 


[m] 


1.9 


0.7 


height 


[m] 


5.0 


1.47 


source 


height [m] 


2.7 2.38 


0.72 



Table 2: Ratios R of measured and predicted ^^Ge production rates in the two GALLEX 
51 Cr radioactive source experiments, Crl [15] and Cr2 |T6], and the SAGE ^^Cr [T71 [18] 
and '^''Ar [13] radioactive source experiments, as reported in Ref. [13]. We give also the 
radii and heights of the GALLEX and SAGE cylindrical detectors and the heights from 
the base of the detectors at which the radioactive sources were placed along the axes of 
the detectors. 

and 37 Ar [TB] radioactive source experiments, as reported in Ref. [13j, are listed in Tab. [21 
Since the weighted average, 

i?Ga = 0.88 ± 0.05 , (2) 

is smaller than unity by more than 2a, it can be interpreted as an indication of the 
disappearance of electron neutrinos due to neutrino oscillations [22l [231 El]- The 
in the absence of oscillation is 8.19 for 4 degrees of freedom, corresponding to a 8.5% 
goodness-of-filEl, as shown in Tab. [31 Therefore, a fluctuation of the data in the case of 
no oscillations cannot be excluded. However, since from a physical point of view it is 
interesting to explore possible indications of non-standard physics, in the following we 
consider the case of neutrino oscillations. 

In the effective framework of two-neutrino oscillations, which is appropriate in the 
case of short-baseline oscillations generated by a squared-mass difference much larger 

data. Similar results have been published recently in a PhD thesis i29j and discussed at the Neutrino 
2008 Conference [30 : i?(Crl) = 0.997 ± 0.11 and i?(Cr2) = OMTto H in a standard rise-time analysis; 
i?(Crl) = 0.953 ± 0.11 and i?(Cr2) 0.812l[J i? in a pulse-shape analysis. We have verified that our 
results are stable against such small changes of the data. 

The goodness-of-fit is the probability to obtain a worse fit under the assumption that the model 
under consideration is correct (see Ref. [3T]). It is the standard statistic used for the estimation of the 
quality of a fit obtained with the least-squares method, assuming the validity of the approximation in 
which Xmin ^ distribution with NDF — iVn — iVp degrees of freedom, where is the number of 
data points and iVp is the number of fitted parameters. The fit is usually considered to be acceptable if 
the goodness-of-fit is larger than about 1%. 
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PG 


NDF 






2 


1 


1 


3 




GoF 






0.23 


0.44 


0.057 


0.28 



Table 3: Values of Xmin^ number of degrees of freedom (NDF) and goodness-of-fit (GoF) 
for the fit of different combinations of the results of the Gallium radioactive source exper- 
iments and the Bugey and Chooz reactor experiments. The first three lines correspond 
to the case of no oscillations (No Osc). The following five lines, including the best-fit 
values of sin^ 2{} and Am^, correspond to the case of oscillations (Osc). The last three 
lines describe the parameter goodness-of-fit (PG) [32j . 



than AmgoL and Am^rpj^ (see Refs. [HUH]), the survival probability of electron neutrinos 
and antineutrinos with energy E,y at a distance L from the source is given bjfl 

UL, E,) = l~ sin^ 2^ sin^ ( , (3) 



V 4^. 

where is the mixing angle and Am^ is the squared-mass difference. The fit of the data 
gives information on the values of the mixing parameters sin^ 2{} and Am^ . 

In our calculation, the theoretical value of the ratio R of the predicted ''^Ge production 
rates in each of the Gallium radioactive source experiments in the cases of presence and 
absence of neutrino oscillations is given by 

E.(B.R.).crJdV^L-2 ' 

where i is the index of the z/g lines emitted in ^^Cr or ^''Ar, which are listed in Tab. [H 
The measured ratios are listed in Tab. [21 together with the dimensions of the detectors, 
which we approximate as cylindrical, and the height from the base of each detector at 
which the radioactive sources were placed along the axis of the respective detector. We 
averaged the neutrino path length L with a Monte Carlo integration over the volume V 
of each cylindrical detector. 

In the separate analysis of the result of each Gallium radioactive source experiment in 
terms of neutrino oscillations, the two mixing parameters cannot be determined through 



The symmetry under CPT transformations, which is a characteristic of aU rclativistic local quan- 
tum field theories, implies that the survival probabilities of neutrinos and antineutrinos are equal (see 
Ref. [11] )■ 
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■ - 99.73% C.L. (3a) 
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■ - 99.73% C.L. (3a) 
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99.73% C.L. (3a) 



Figure 1: Allowed regions in the sin^ 2'd-/S.'w? plane obtained from the fits of the results of 
the two GALLEX ^^Cr radioactive source experiments, Crl and Cr2, and the SAGE ^^Cr 
and ^''Ar radioactive source experiments. The curves in the GALLEX Crl and SAGE 
^^Cr plots exclude the region on the right. In the GALLEX Cr2 and SAGE ^^Ar plots, 
the pairs of la and 2a curves delimit allowed regions, whereas the 3cr curves exclude the 
region on the right. 



a least-squares analysis from one data point. Therefore, we adopt a Bayesian approach, 
as done in Ref . [33] , considering i? as a random variable with a uniform prior probability 
distribution between zero and one. Then, if i?obs is the observed value of i?, the normalized 



posterior probability distribution of R is given by 

^ ' ^ j,'dRp{R^^,\R) ^ ' 

Here, p{Rohs\R) is the sampling distribution of -Robs given R, which we assume to be 
Gaussian with standard deviation equal to the experimental uncertainty. The allowed 
interval of R with a given Bayesian Confidence Level is given by the Highest Posterior 
Density interval with integrated probability equal to the Confidence Level. Figure [T] 
shows the resulting allowed regions in the sin^ 2'i9-Am^ plane. One can see that the first 
GALLEX source experiment (Crl) and the ^^Cr SAGE source experiment, in which the 
measured rate is within la from unity, imply only upper limits for the mixing parameters. 
On the other hand, the analyses of the second GALLEX source experiment (Cr2) and 
the ^''Ar SAGE source experiment give 2(T allowed bands, which have a large overlap for 
Am2 > leVl 

Let us now discuss the combined fit of the four Gallium source experiments. Since 
there are enough data points to determine the two mixing parameters sin^ and Am^, 
we abandon the Bayesian approach in favor of a standard frequentist least-squares fit. 
This method is based on a global minimization of the iii the sin^ 2d-^m? plane and 
the calculation of the Confidence Level contours corresponding to a Ax^ with two degrees 
of freedom: Ax^ = 2.30,6.18,11.83 for 68.27% (la), 95.45% {2a) and 99.73% (3cr) C.L., 
respectively (see Ref. [31]). 

The result of the combined least-squares analysis of the four Gallium source experi- 
ments is shown in Fig. [21 One can see that there is an allowed region in the sin^ 2'd-/S.rn? 
plane at la for Am^ > 0.6 eV^ and 0.08 < sin^ 2'(9 < 0.4. The values of Xmin^ the number 
of degrees of freedom (NDF), the goodness-of-fit (GoF) and the best-fit values of the 
mixing parameters are given in Tab. [3l The value of the the goodness-of-fit (23%) shows 
that the fit is acceptable. 

Table m shows the allowed ranges of sin^ 2t? and Am^ obtained from the corresponding 
marginal Ax^ = — x^m ^ig. O The presence of 2a lower limits for sin^ 2'd and Am^ 
in spite of the absence of a 2a lower limit in the sin^ 2'd-/S.rn? plane in Fig. [2] is an effect 
due to the statistical analysis: for one parameter 2a corresponds to Ax^ = 4, whereas 
for two parameters it corresponds to Ax^ = 6.18. Hence, it is fair to conclude that there 
is an indication of a possible neutrino disappearance due to neutrino oscillations with 
s\.v?2'd > 0.03 and /S.rn? > 0.1 eV^ at a confidence level between one and two sigmas 
(- 70 - 90%C.L.). 



3 Bugey 

The disappearance of electron antineutrinos have been investigated by several reactor 
neutrino experiments at different baselines (see Refs. [3ll[TT]). Since, according to Eq. ([3]), 
the survival probabilities of neutrinos and antineutrinos are equal, the interpretation 
of the results of Gallium radioactive source experiments in terms of electron neutrino 
disappearance can be compared directly with the results of reactor neutrino experiments. 

In this section we consider the results of the reactor short-baseline Bugey experiment 
[25] , which put the most stringent constraints on the disappearance of electron antineu- 
trinos due to Am^ > O.leV^. 
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Table 4: Allowed ranges of sin^ 2{} and Am^ from the combined fit of the results of 
Gallium radioactive source experiments, from the fit of the results of the Bugey reactor 
experiment, and from the combined fit. The dash indicates the absence of limits. 



Reactor neutrino experiments detect electron antineutrinos through the inverse neu- 
tron decay process 

i^e + P ^ n + . (6) 
The neutrino energy and the positron kinetic energy Tg are related by 

E^ = Te + Tr, + me + mn-mpC:^Te + 1.8 MeV , (7) 

where is the negligibly small recoil kinetic energy of the neutron. In the Bugey ex- 
periment the survival probability of electron antineutrinos was measured at three source- 
detector distances: Lj = 15,40,95 m, for j = 1,2,3, respectively. We use the ratio of 
observed and expected (in the case of no oscillation) positron spectra given in Fig. 17 of 
Ref. p5|, in which there are Nj = 25,25, 10 energy bins. We analyze the data with the 
following x^, taken from Ref. [25] : 



,2 

2 ' 



where Eji is the central energy of the ith bin in the positron kinetic energy spectrum 
measured at the Lj source-detector distance, i?^^^ and i?*^^ are, respectively, the corre- 
sponding measured and calculated ratios. The uncertainties aji include the statistical 
uncertainty of each bin and a 1% systematic uncertainty added in quadrature, which 
takes into account the uncertainty of the spectrum calculation (with a total of about 5% 
uncorrelated systematic uncertainty over 25 bins). The coefficients {Aaj + b{Eji — Eq)), 
with Eq = 1 MeV, were introduced in Ref. [25] in order to take into account the sys- 
tematic uncertainty of the positron energy calibration. The value of as a function of 
sin^ 2{} and Am^ is calculated by minimizing Eq. ([H]) with respect to the five parameters 
A, ttj (j = 1,2,3), 6, which have, respectively, uncertainties aA = 0.048, cTq^ = 0.014 
(Tb = 0.02 MeV~^ [25]. Following Ref. [35], we approximate the neutrino fiux, the detec- 
tion cross section and the detection efficiency as constants in each energy bin. Then, i?*-^*^ 
is given by 

^ / IS-^E^I: rr ^T. F{E, T,) P,^^,^ m E^) 
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Figure 2: Allowed regions in the sin^ plane and marginal Ax^'s for sin^ 2'd and 

Am? obtained from the combined fit of the results of the two GALLEX ^^Cr radioactive 
source experiments and the SAGE ^^Cr and ^^Ar radioactive source experiments. The 
best-fit point corresponding to Xmin is indicated by a cross. 



Here Tg and E^, are, respectively, the positron kinetic energy and the neutrino energy, 
related by Eq. ([7]), whereas E is the measured positron kinetic energy, which is connected 
to Te by the energy resolution function of the detector F{E, Tg). We considered a Gaussian 
energy resolution function with standard deviation 0.252 A/-E/4.2MeV MeV [25]. The 
quantities /S.Ej are the widths of the energy bins in each detector. The integration over 
the neutrino path length L is performed by a Monte Carlo which takes into account the 
geometries of the reactor and of the detectors and their relative positions [36j . 
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Figure 3: 90% C.L. exclusion curves in the sin^ 2t9-Am^ plane obtained from a raster-scan 
analysis of Bugey data (solid line) and from a standard global least-squares fit (dashed 
line) . 



With this method we obtained the 90% C.L. raster-scaij^l exclusion curve shown in 
Fig. [21 which is similar to the original 90% C.L. raster-scan Bugey exclusion curve in 
Ref. [25]. 

Let us emphasize that the raster-scan method is statistically weak, because it does 

^ In the raster-scan method, Xmin found for each fixed value of Am^. The corresponding upper 
hmit for sin^ 2i9 is calculated as the value of sin^ 2i? for which the cumulative distribution function of 
Ax^ = ~ Xmim which has one degree of freedom, is equal to the Confidence Level (Ax^ = 2.71 for 
90% C.L.). 



9 



00 -■ 



CD -■ 

CM 




Figure 4: Allowed regions in the sin^2'(9-Am^ plane and marginal Ax^'s for sin^ 2'(9 
and Am^ obtained from the least-squares analysis of Bugey data. The best-fit point 
corresponding to xLin is indicated by a cross. 



not have proper coverage [37] . We presented in Fig. [3] the raster-scan exclusion curve 
only to show by comparison with the analogous figure in Ref. p5] that our analysis of the 
Bugey data is acceptable. The dashed line in Fig. [H] shows the 90% C.L. Bugey exclusion 
curve obtained with the standard least-squares method, which we adopted also in the 
previous Fig. [2] and the following Figs. HHHl From Fig. [3] one can see that the 90% C.L. 
raster-scan exclusion curve overcovers for all values of Am^, except for small intervals 
around Am? ~ 0.9 eV^ and Am^ ~ 1.9 eV^. 

Figure m shows the allowed regions in the sin^ plane and the marginal Ax^'s 
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Figure 5: Best fit of Bugey data (points witli error bars [2S])- Tlie tliree panels sliow 
tlie ratio R of observed and expected (in tlie case of no oscillation) event rates at the 
three source-detector distances in the Bugey experiment as functions of the measured 
positron kinetic energy E (see Eq. (jO])). In each panel, the solid and dashed histograms 
correspond, respectively, to the best-fit values of {Aaj + b {Eji — Eq)) E^^^ and B}jf (see 
Eq. dHD). 



for sin^ 2'd and Am? obtained from the least-squares analysis of Bugey data. The value 
and location in the sin^ 2'd- Aw? plane of the minimum of the x^? the number of degrees of 
freedom (NDF) and the goodness-of-fit (GoF) are given in Tab. |3l The fit is satisfactory, 
since the goodness-of-fit is 67%. The best-fit value of the oscillation parameters and the 
small la allowed regions in Fig. H] are in favor of neutrino oscillations. However, the 2(T 
and 3cT contours in Fig. H] provide only upper limits to neutrino oscillations. Also, the 
value of the iii the case of absence of oscillations and the corresponding goodness-of-fit 
(63%) do not allow us to exclude the absence of oscillations. 

The reason of the hint in favor of neutrino oscillations given by the Bugey data is 
illustrated in Fig. El where the histogram relative to the best fit is shown against the 
Bugey -R^P's. With the help of the histogram, one can see that there is a weak hint 
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Figure 6: Allowed regions in the sin^ 2{}-Am'^ plane and marginal Ax^'s for sin^ 2{} and 
Am^ obtained from the combined fit of the results of the two GALLEX ^^Cr radioactive 
source experiments, the SAGE ^^Cr and '^''Ar radioactive source experiments and the 
Bugey reactor experiment. The best-fit point corresponding to Xmin is indicated by a 
cross. 



of oscillations. The la allowed regions in Fig. H] have very narrow Am^ ranges around 
0.9 eV^ 1.85 eV^ and 3eV^ because slight shifts of Am^ from these optimal values spoil 
the agreement with the data of the histogram in Fig. O 

Table H] shows the marginal allowed ranges of sin^ 2'd and Am^ obtained from the 
corresponding A^^'s in Fig. |H One can see that there is a hint of neutrino oscillations 
with 0.02 < sin^ 2'd < 0.08 and Am^ ^ 1.8 eVl 
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From a comparison of Figs. [2]and|l]one can see that the allowed regions of the Gallium 
radioactive source experiments and the Bugey experiment are marginally compatible for 
sm 2d ~ 0.1 and Am^ > leV^ Fi gure [6] shows the allowed regions obtained from the 
combined fit. Since the Bugey data are statistically dominant, the curves in Fig. [H] are 
not very different from those in Fig. HI which have been obtained from the fit of the 
Bugey data alone. The inclusion of the Gallium data has the effect of eliminating the 
la allowed region at Am^ 0.9 eV^ and of disfavoring at Icr values of siv? 2^} smaller 
than about 2 x 10^^. The value and location of Xmin; ^^e number of degrees of freedom 
and the goodness-of-fit are listed in Tab. [31 One can see that the Gallium data do not 
spoil the good fit of the Bugey data. Indeed, the value of the parameter goodness-of-fit@ 
[H2] reported in Tab. [H] shows that the Bugey and Gallium data are compatible under 
the hypothesis of neutrino oscillations. The marginal allowed ranges of sin^ 2^9 and Am^ 
obtained from the corresponding Ax^'s in Fig. [6] are given in Tab. HI 



4 Chooz 

In this section we consider the result of the long-baseline reactor neutrino experiment 
Chooz |26j, which gives limits on neutrino oscillations which are comparable with those 
of the Bugey experiment for Am^ > 2eV^. 

In the Chooz experiment the ratio of the number of observed events and that expected 
in the absence of neutrino oscillations is 

i?chooz = 1.01 ±0.04. (10) 

The value of this ratio puts a constraint on the disappearance of electron (anti)neutrinos 
with energies in the MeV range at distances smaller than about 1 km. This corresponds 
to a constraint on sin 2§ for Am2 > lO'^eVl In the range of sensitivity of the Gallium 
radioactive source experiments, Am^ > 10^^ eV^ (see Figs. [2]), the oscillation length of 
reactor antineutrinos is much shorter than the Chooz source-detector distance. In this 
case, the Chooz experiment is only sensitive to the averaged survival probability 

{Pi-) (-)) = 1 - ^ sin^ 2^ . (11) 

Therefore, the Chooz result in Eq. ffTUj) can be combineclll with the results of the Gallium 
radioactive source experiments simply by considering it as a measurement of sin^ 2'd: in 
the Bayesian approach of Eq. ([5]) 

sin^ 2^9 < 0.071, 0.15, 0.23, (12) 

at 68.27% (Id), 95.45% {2a), 99.73% (3(t) Bayesian Confidence Level, respectively. 

^ The value of (Ax^;jj)a+b corresponding to the parameter goodness-of-fit of two experiments A and 
B is given by (Xmin)A+B — [(Xmin)A + (Xmin)B] ■ It has a distribution with number of degrees of freedom 
NDF = Pa + -Pb — -Pa+b, where Pa, -Fb and Pa+b are, respectively, the number of parameters in the 
fits of A, B and A-l-B data [52] , 

^ In our figures we considered Ato^ in the range 10~^ — 10^ eV^. For simplicity, we neglected the 
small Am^ dependence of the CHOOZ exclusion curve for Am2 < 4 X IQ-^eV^ (see Fig. 55 of Ref. 
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Table 5: Allowed ranges of sin^ 2'd and Am^ from the combined fit of the results of 
the Bugey and Chooz reactor experiments, the Gallium radioactive source and Chooz 
reactor experiments, and the Gallium radioactive source and the Bugey and Chooz reactor 
experiments. The dash indicates the absence of limits. 



First, we performed a combined frequentist least-squares analysis of the Bugey and 
Chooz data, which yielded the allowed regions in the sin^ 2t9-Am^ plane shown in Fig. [71 
the best fit values of the mixing parameters reported in Tab. [3l and the marginal allowed 
ranges listed in Tab. [5l One can see that the addition of the Chooz result to the Bugey 
data analysis has the effect of improving slightly the upper limit on sin^ 2'd for Am^ > 
3 eV^ and that of excluding values of sin^ 2'd larger than about 0.1 for Am^ ^ 3 x 10~^ eV^, 
where Bugey is not sensitive. In the intermediate range of Am^, where Bugey is sensitive 
to the oscillations, the addition of the Chooz result weakens the hint in favor of oscillations 
given by the Bugey data: the la allowed regions in Fig. H] are stretched towards small 
values of sin^ 2{} in Fig. [71 However, the best-fit value of the mixing parameters remain 
unchanged, because of the dominance of the Bugey data. From Tab. [SI one can see that 
the parameter goodness-of-fit implies that Bugey and Chooz results are compatible under 
the hypothesis of neutrino oscillations, but the goodness-of-fit obtained in the case of no 
oscillations do not allow us to exclude this possibility. 

From the comparison of Eq. ( [T2l) and Fig. [21 one can see that the results of the Chooz 
and the Gallium radioactive source experiments are compatible only at the 2a level. In 
fact the parameter goodness-of-fit reported in Tab. [3 shows a tension between Gallium 
and Chooz data under the hypothesis of neutrino oscillations. 

Figure M shows the allowed regions in the sin^ 2t9-Am^ plane obtained with the com- 
bined least-squares fit of Gallium and Chooz data. The values of Xmm Goodness 
of Fit and the best-fit values of the mixing parameters are given in Tab. [31 It is clear 
that the combined fit is not good, since the results of Chooz and the Gallium radioactive 
source experiments are in contradiction regarding neutrino disappearance. The marginal 
allowed ranges of sin^ 2{} and Am^ in Tab. [5lare of little interest, since the minima of the 
corresponding Ax^'s in Fig. [HI are very shallow, except for the upper bound on sin^2'i9 
driven by Chooz data. As one can see from the allowed regions in the sin^ 2'd-Am'^ plane 
in Fig. [HI the Chooz bound on sin^ 2'i9 in Eq. fll2l) is weakened by the results of the Gal- 
lium radioactive source experiments in a significant way only for Am^ > 10~^ eV^ at the 
la level. 

Finally, we performed a combined fit of the results of Bugey, Chooz, and Gallium data. 
The resulting allowed regions in the sin^ 2'(9-Am^ plane are shown in Fig. [9l The best fit 
values and the marginal allowed ranges of the mixing parameters are listed, respectively. 
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Figure 7: Allowed regions in the sin^ plane and marginal Ax^ 's for sin^ 

and Am^ obtained from the combined fit of the results of the Bugey and Chooz reactor 
experiments. The best-fit point corresponding to Xmin is indicated by a cross. 



in Tabs. [3] and [5l One can see that the Gallium and Chooz data tend to compensate each 
other, leading to results which are similar to those obtained in the analysis of Bugey data 
alone. The value of the parameter goodness-of-fit reported in Tab. [S] does not allow us to 
exclude the compatibility of the Bugey, Chooz and Gallium data under the hj^othesis of 
neutrino oscillations. Also the goodness-of-fit obtained in the case of no oscillations, given 
in Tab. [3], is acceptable. Therefore, we can only conclude that the combined analysis of 
all the experimental data that we have considered is compatible both with the case of no 
oscillations and with the hint in favor of neutrino oscillations with 0.02 < sin^ 2'd < 0.08 
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10"^ 10"^ 10"' 2 4 6 8 10 



Figure 8: Allowed regions in the sin^ 2'd-A.w? plane and marginal Ax^'s for sin^ 2'd and 
Am^ obtained from the combined fit of the results of the two GALLEX ^^Cr radioactive 
source experiments, the SAGE ^^Cr and ^^Ar radioactive source experiments and the 
Chooz reactor experiment. The best-fit point corresponding to x^^-^ is indicated by a 



and Am^ 1.8 eV^ found in the analysis of Bugey data. 

5 Conclusions 

We interpreted the deficit observed in the Gallium radioactive source experiments as a 
possible indication of the disappearance of electron neutrinos. We have analyzed the data 
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Figure 9: Allowed regions in the sin^ 2i?-Am^ plane and marginal Ax^'s for sin^ 2f? and 
Am^ obtained from the combined fit of the results of the two GALLEX ^^Cr radioac- 
tive source experiments, the SAGE ^^Cr and ^^Ar radioactive source experiments and 
the Bugcy and Chooz reactor experiments. The best-fit point corresponding to xLin 
indicated by a cross. 



in the effective framework of two-neutrino mixing, which describes neutrino oscillations 
due to a Am^ that is much larger than the solar and atmospheric ones. We found that 

there is an indication of electron neutrino disappearance due to neutrino oscillations with 
sin^ 2^ > 0.03 and Am^ ^0.1 eV^. We have also studied the compatibility of the data of 
the Gallium radioactive source experiments with the data of the Bugey and Chooz reactor 
short-baseline antineutrino disappearance experiments in the same effective framework of 
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two-neutrino mixing, in which the disappearance of neutrinos and antineutrinos are equal. 
We found that the Bugey data present a hint of neutrino oscillations with 0.02 < sin^ 2'd < 
0.08 and Am^ «i 1.8 eV^, which is compatible with the region of the mixing parameters 
allowed by the analysis of the data of the Gallium radioactive source experiments. We 
have also performed combined analyses of the Bugey and Chooz data, of the Gallium and 
Bugey data, of the Gallium and Chooz data, which show that the Bugey and Chooz data 
are compatible, the Gallium and Bugey data are compatible, and the Gallium and Chooz 
data are marginally compatible. The weak indication in favor of neutrino oscillations 
found in the analysis of the Bugey data persists in the combined analyses of the Bugey 
data with the Gallium and Chooz data. However, we cannot exclude the absence of 
oscillations. 

From a physical point of view, a hint in favor of short-baseline neutrino oscillations 
generated by Am^ > 0.1 eV^ is extremely interesting. This squared-mass difference is 
too large to be compatible with the three-neutrino mixing scheme inferred from the 
observation of neutrino oscillations in solar, very-long-baseline reactor, atmospheric and 
long-baseline accelerator experiments, in which there are only two independent squared- 
mass differences, AmgQL ~ 8 x 10~^eV^ and Am^rpj^ 3 x 10"'^ eV^. Therefore, the 
results of our analysis indicate the possible existence of at least one light sterile neutrino 
Vg (see Refs. O [6l [H]). We think that it is very important to explore this intriguing hint 
of new physics beyond the Standard Model. 

As already discussed in Ref. [21], short-baseline transitions have an influence 

on the interpretation of all experiments with an initial beam. In the existing solar and 

atmospheric neutrino experiments the survival probability of is the averaged one in 
Eq. ffTTl) . However, the uncertainties of the experimental data and our knowledge of the 

initial flux do not allow us to exclude — > transitions at the level of about 20% in 
the case of solar neutrinos [38j and about 30% (see Ref. |llj) in the case of atmospheric 
neutrinos. 

Future experiments which are well suited for finding small transitions are 

those with a source producing a'l/g flux which is known with high accuracy. Since sterile 
neutrinos are invisible, — transitions can be revealed either by measuring a disap- 
pearance of Ve's without or V production in the detection process or by measuring a 
disappearance of'z/g's due to oscillations with a squared-mass difference much larger than 
AmgQL ^^ATM- aware of the following possibilities: Beta-Beam experiments 

[39] which have a pure z/g or z/g beam from nuclear decay (see the reviews in Refs. [lOll^): 
Neutrino Factory experiments in which the beam is composed of Ve and z/^, from /i+ de- 
cay, or z/g and z/^, from /i~ decay (see the review in Ref. [l2l HQ]); Mossbauer neutrino 
experiments, with a Ve beam produced in recoiless nuclear decay and detected in recoiless 
nuclear antineutrino capture [13]; the LENS detector [HHIH] with an artificial Megacurie 
z/g source [10]. Let us also notice the very interesting possibility to reveal the existence 
of sterile neutrinos in the flux of high-energy astrophysical neutrinos after their passage 
through the Earth by measuring the peculiar matter effects [171 H8]. 
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